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ABSTRACT 

The statistical mechanics of isolated collisionless self-gravitating systems is a long-held puz- 
zle, which has not been successfully established for nearly 50 yr. In this work, we employ a 
phenomenological entropy form of ideal gas, first proposed by White & Narayan, to revisit 
this issue. By calculating the first-order variation of the entropy, subject to the usual mass- 
and energy-conservation constraints, we obtain an entropy stationary equation. Incorporated 
with the Jeans equation, and by specifying some functional form for the anisotropy parame- 
ter /?, we numerically solve the two equations, and demonstrate that the velocity anisotropy 
parameter plays an important role to attain a density profile that is finite in mass, energy, and 
spatial extent. If incorporated again with some empirical density profile from simulations, our 
theoretical predictions of the anisotropy parameter, and the radial pseudo-phase-space den- 
sity p/a^ in the outer non-gravitationally degenerate region of the dark matter halo, agree 
rather well with the simulation data, and the predictions are also acceptable in the middle 
weak-degenerate region of the dark halo. The disagreements occur just in the inner strong- 
degenerate region because of the ignoring of gravitational degeneracy. As far as we know, our 
results may be the first theoretical predictions based on the entropy principle that can partially 
match the empirical data. The second-order variational calculus reveals the seemingly para- 
doxical but actually complementary consequence that the equilibrium state of self-gravitating 
systems is the global minimum entropy state for the whole system under long-range violent 
relaxation, but simultaneously the local maximum entropy state for every and any small part 
of the system under short-range two-body relaxation and Landau damping. This minimum- 
maximum entropy duality means that the standard second law of thermodynamics needs to 
be re-expressed or generalized for self-gravitating systems. We believe that our findings, es- 
pecially the complementary second law of thermodynamics, may provide crucial clues to the 
development of the statistical physics of self-gravitating systems as well as other long-range 
interaction systems. 

Key words: methods: analytical - cosmology: theory — dark matter — large-scale structure 
of Universe. 



1 INTRODUCTION 

Cosmological TV-body simulations of cold dark matter have 
revealed many empirical relationships concerning collisionless 
self-gravitating dark matter haloes, such as the density profile, 
the pseudo-phase-space density profile, the distributions of 
shape and spin parameters, the velocity- ani sotr opy-parameter 
profil e and the distribution o f velocity dispersion jNavarro et al.l 
I2OO8I : IWang & Whit3 |2009|) . They used to be considered as 
being universal, i.e. different cold dark matter haloes obey 
the same empirical laws only with different scale parameters, 
which are independent of the halo mass, the initial density 
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fluctuation spectrum and the values of the cosmological param- 
eters. The famous one of these universalities is the universal 
density profile of cold dark matter halo es. The se universal pro- 
files are usually e xpressed as NFW jNavarro, Frenk & White 



1 19951 1 19961 1 19971: iBuUock et alj I2OO ih or iMoore et all l\ 999 
or even other specific forms, such as the lEinastd ( 119651) pro- 
file. We should notice that, however, more and more recent 
high-resolution A^'-body simulations claim that the dark halo 
density profiles are not strictly universal, which may depend 
on ha l o mass, merger history or substruct ur es (see Ijing & Sutd 



20051: IM erritt et al. 
iNavarro et al.,i2008i) . 
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In fact, long before these simulation results of dark mat- 
ter haloes were found, astronomical observations had already re- 
vealed that some self-bounded gravitating systems exhibit remark- 
able structural regularities. For example, elliptical galaxies display 
a quas i-universal luminosity profile described bv Ide VaucouleursI 
( Il948h 's R^^*^ law, and gl obular clusters are generall y well fitted 
by the Michie-King model teinnev & Tremaingl2008l) . 

Although these empirical relationships are simple and prac- 
tical, we just know what, but not why they behave in such 
ways, since they have no physical basis. Up to now, however, 
few of the physical origins concerning these empirical results 
have been investigated, except for the universal dark halo den- 
sity profiles, whose orig in was explo red as soon as the universal- 
ity was found tSver & White! [l993) . Thereafter, many other ex- 
planations for the origin of dark halo density pro files were pro- 
posed , and most of them, similar to the one bv lSver&Whia 
( Il998h . are concentrated on various dynamical processes, such 
as collapsing, mergers of subhaloes or accretion of surround 
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Ascasibaretal.'2004; Williams, Babul & Dalc antonl2004lLu et al.l 
2006 ; Bellovary et al. 2008 ; El-Zant 2008). However, a recent work 
by IWang & White! (2009) suspected this kind of dynamical ori- 
gins. According to their simulations of hot dark matter haloes, they 
claimed that mergers do not play a pivotal role in establishing the 
universalities, and thus their results contradict the models which 
explain the universalities as consequences of mergers. 

Is it possible that all these empirical relationships originate 
from statistical mechanics, or more concretely, do they reflect 
the properties of self-gravitating systems in statistical equilibrium 
states? If so, from theories of thermodynamics and statistical me- 
chanics, we know that an equilibrium state of a thermodynamic sys- 
tem does not depend on any specific dynamical processes or mech- 
anisms of how to attain it. If dark matter haloes really attain the 
statistical equilibrium states, then it would be very easy to under- 
stand the finding bv I Wang & White! ( i200^ . why the universalities 
are irrelevant to mergers. 

The investigation of equilibrium statistical mechanics of iso- 
lated collisionless self-gravitating syste ms is a lo ng-held issue, 
which can date back to the early 1960s. lAntonov! (1962) proved 
that global entropy maxima do not exist for self-gravitating sys- 
tems, which might be the first study that exhibits the great com- 
plexity of_the_statistic^_mechanics for self-gravitating systems. 
Later, !Lvnden-Bell! jl967 ) in his seminal work proposed 'vio- 
lent relaxation' to account for the formation process of ellipti- 
cal galaxies. Just in the same work, he derived a 'most proba- 
ble' distribution of a system in its 'equilibrium state' that has in- 
finite mass, energy and sp atial extent, and hen ce is not reali s tic an d 
not acceptable. F ollowing Lvnden-Bell ('l96y)'s work, Shu (1978), 
!Ku11, Treumann & Bohringer (1997) and Nakamura (2000) , with 
increasing complexity of phase-space partitions, reworked for their 
own coarse-grained distribution schemes. These authors aimed to 
remedy the flaws in the previous works, the most significant of 
which is the mass segregation problem, i.e. in a multi-species com- 
position, the heaviest objects tend to concentrate more centrally 
whereas the lightest ones tend to go towards the border. However, 
these authors did not overcome the problem of infinite mass that 
arises from Lynden-Bell ( 196 7)'s s tatistical approach. We refer the 
reader to IBindoni & Seccol ( 120081) for the comprehensive review 
about the history of these works. 

Since the concept of entrop y is the basis of statistical mechan- 
ics, ever since !Antonovl jl962t) . who used the usual Boltzmann- 



Gibbs entropy, various alternative entropy forms have been at- 
tempted to account for the stability and equilibrium configu- 
rations of self- gravitating systems, s uch as the generalized H- 
function ( Tremaine. Henon & Lvnden- Bell 1986), t he entropy of 
ideal gas ([w hite & N aravan!!l987!) or Tsallis entropy jTsalhsll98^: 
!Hansen etalj ,2005) . There are also other authors using entropy 
arguments to draw conclusions about the final states of vio- 
lent r elaxation (e.g. !Stiavelli & B ertin 19871: !Spergel & Hemauis3 
All these entropies are maximized, subject to the 



usual mass- and energy-conservation constraints |Antonov 19621 : 
!Lvnden-Belll!l967! : !Tremaine etZl!l986! : !white & Naravan!ll987!) . 
or other additional con ditions, such as hydrostatic equilibrium 
( IWhite&Naravan!!l987h . 

Usually, the equilibrium state of a thermodynamic system is 
defined by employing the principle of maximum entropy, which 
states as follows: the realistic thermodynamic relations of the sta- 
ble equilibrium state for any thermodynamic system may be de- 
rived by seeking the probability density function that maximizes 
the Boltzmann-Gibbs entropy (or other plausible entropy forms), 
subject to all relevant constraints. However, the principle of max- 
imum entropy for self-gravitating systems has been continually 
questioned for the past several decades. For example, there is a 
famous argument, originated by J. Binney, against the feasibility 
of Boltzmann-Gibbs entro py, in that such an entr opy form is un- 
bounded above ([s innev & Tremaine!!2008! : see also Tremaine et al.! 
!l986l:!White & Naravarj|19^ , which resembles lA"ntonov| jl962Fs 
earlier result that there exists no global maximum entropy for self - 
gravitating sy stems. Confirming an earlier study bv !Sricfcail ( fl987h , 
!Soker! ( ! 19961) showed that the general iif -functions during the pro- 
cess of violent relaxation, may not necessarily increase at all times 
and may decrease if the system has specific, far from equilib- 
rium, conditions. They further proposed to distinguish between 
two phases of violent relaxation, i.e. in one phase, the //-function 
may decrease with time, whereas in the other it is always a non- 
decreasing function of time. This result suggests that whether H- 
function is a monotonically increasing function o f time during th e 
violent relaxation process is still an open question ( !Kandrud!l987h . 
or equivalently, the equilibrium state of self-gravitating systems 
should not be simply considered as the maximum entropy state 

(Chavanis 2002).^ 

Worse still, !Arad & Lvnden-Bell! ( !2005h revealed that there 
exists an inconsis tency in theories of violent relaxation with 
!Lvnden-Belll ( !l967h and !Nakamur3 ( !2000!) 's distributions. The in- 
consistency arises from the non-transitivity of these theories, i.e. 
a system that undergoes a violent relaxation, then followed by an- 
other, with an addition of energy for the latter, is not equivalent to 
the one that undergoes the whole process directly from the initial 
to the final state. The non-transitivity is destructive, since the 'equi- 
librium state' is dependent on the evolutionary path of the system, 
which implies that a complete methodology of statistical mechan- 
ics is probably impossible for self-gravitating systems, and should 
be replaced by some dynamical approaches. 

Therefore, after unsuccessful attempts for nearly 50 yr, the 
concept of entropy as well as the methodology of statistical me- 
chanics might have been widely regarded as inviable for self- 
gravitating systems. As a result, we are in such an embarrassing 
situation that we have few applicable analytic theories concerning 
self-gravitating systems, such as elliptical galaxies and dark matter 
haloes, except for the Poisson-Vlasov equations (or the Jeans equa- 
tion) and the fitting formulae from either TV-body simulations or 
observations. 

At the same time, these problems and difficulties agonize 
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statistical physicists as well. They raised such questions for the 
general long-range interaction systems, as non-additivity or non- 
extension, ensemble inequivalence and negative sp ecific heat in 
micro canonical ensemble and ergodicity breaking jCampa et al.l 
Moreover, since it is not known how to define the equilib- 
rium state of a self-gravitating system, such a state that an equili- 
brated system attains is vaguely called the quasi-stationary state. 

It seems that, in order to establish a successful statistical me- 
chanics of self-gravitating systems, some new ideas or even con- 
ceptual revolutions may be desired to get out of this long-term im- 
passe. In this work, we explore whether the puzzle concerning the 
equilibrium statistical mechanics of self-gravitating systems can be 
solved by introducing some new thoughts, or by revising some fun- 
damental laws or principles of physics. 

Our paper is arranged as follows. In Section|2] we firstly pro- 
pose a mean-field model to separate gravitational forces into long- 
and short-range parts, and subsequently identify three distinct re- 
laxation mechanisms, i.e. violent relaxation, two-body relaxation 
and gravitational Landau dampi ng. Then by calcu l ating the first- 
order variation of the entropy of IWhite & NaravanI l ll987ll . we de- 
rive an entropy stationary equation. In Section |3] combined with 
the Jeans equation, we solve the two equations in two cases. First, 
we show that the velocity anisotropy parameter is of great signifi- 
cance to make the density profiles of dark matter haloes have finite 
mass, energy and extent. Then, ignoring the gravitational degen- 
eracy in the inner region of the dark matter halo, we find that the 
predicted anisotropy parameter and the pse udo-phase-space den- 
sity agree rather well with the simulations o f iNavarro et ai] ( I2OO8I) 
in the outer non-degenerate region, and are also in rough agree- 
ments with the simulations in the middle weak-degenerate region 
of the dark halo. An upturn deviation from the power law of the 
pseudo-phase-space density at the ha lo outskirts is also p redicted, 
which is consistent with the finding bv lLudlow et al. In Sec- 

tion |4l we generalize the standard second law of thermodynamics 
to long-range self-gravitating systems, which states in the comple- 
mentary way as follows: for an isolated self-gravitating system, the 
entropy can automatically decrease under long-range violent relax- 
ation, tmtil violent relaxation terminates and the system attains the 
equilibrium state with the global minimum entropy for the whole 
system; yet, in this very equilibrium state, the system is simulta- 
neously in the local maximum entropy state for every and any part 
of the system under short-range two-body relaxation and Landau 
damping. In Section |5] we conclude that, with this revised, com- 
plementary second law of thermodynamics, the concept of entropy 
is still valid, and thus in principle, the equilibrium statistical me- 
chanics can also be established for such long-range self-gravitating 
systems. 



2 BASIC THEORIES AND FORMULAE 

2.1 Mean-fleld model for self-gravitating systems 

It is us ually believed that, undergoing the violent relaxation pro- 
cedure jLvnden-Belill 19671) , a s elf-gravitating syste m can attain a 
so-called quasi-stationary state dCampa et alJlaOOgh . i.e. the virial 
equilibrium state. Throughout this paper, we tentatively call this 
quasi-stationary state as equilibrium state. 

The appropriate gravitational model should be the mean-field 
model, which is the first step towards the correct statistical me- 
chanics of self-gravitating systems. The following is a brief intro- 
duction to the idea of mean-field model, which was suggested by 
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Figure 1. (a) Illustration of the separation of gravitational force. The system 
possesses some geometrical symmetry in the equilibrium state, so that the 
long-range parts of gravitational forces of the particles can be collectively 
manifested as a mean-field force, with short-range residual effects deter- 
mined by nearby particles, (b) The three relaxation mechanisms of self- 
gravitating systems, i.e. violent relaxation, two-body relaxation and Landau 
damping. The system is driven by violent relaxation into the virial equilib- 
rium state, which is maintained by Landau damping and two-body relax- 
ation. 



ISvlos Labinl ( |2008|) . For the detailed discussion of the mean-field 
model for self-gravitating systems, we refer the reader to their orig- 
inal article. There are also many discussions about the mean-field 
model that can be found in Binney & Tremaine (2008). 

Usually, the standard statistical mechanics deals with short- 
range interaction systems. The typical feature of these systems is 
that the interaction between two particles is strong and repulsive 
when they are very close to each other, but on large distances no 
force is exerted between two particles. As a result, particles most 
of time move at nearly constant velocity, and then are subject to 
violent and instantaneous accelerations when they collide with one 
another. In this situation, the typical relaxation mechanism arises 
from two-body collisions. 

One of the difficulties of developing the statistical mechan- 
ics of self-gravitating systems lies in that gravitational interaction 
is long-range attractive force, which means that a particle is cou- 
pled to almost all other particles at all spatial scales of the system. 
So it seems that all different range of scales are equally impor- 
tant to the dynamical properties of self-gravitating systems. For- 
tunately, we can circumvent this obstacle by constructing a gravi- 
tational model, in which the long-range part of gravitational force 
between the particles is collectively manifested as a mean field for a 
finite self-gravitating system, and the short-range part arising from 
the nearby particles is treated as fluctuations in this mean gravita- 
tional field. Note that the equilibrated self-gravitating systems usu- 
ally possess some geometrical symmetry, by virtue of which the 



4 P. He and D.-B. Kang 



long-range coupling between the particles can be broken when the 
system is in its equilibrium states. 

We illustrate this separation of gravitational force with a sim- 
plified one-dimensional model as follows. See Fig. [TJa), for the 
one-dimensional system in equilibrium, the particles are homoge- 
neously distributed in the mean field, indicated by the dashed line, 
with a mean separation as Is - In this model, if a particle is initially 
located at point A, it will not feel any gravitation. When the par- 
ticle moves towards the right side, however, it will be more and 
more affected by the gravitational force of particle Pi. When the 
particle crosses over Pi on the right side, continuously moving far- 
ther and farther away, the gravitational force exerted on this par- 
ticle becomes less and less with the increasing distance from Pi, 
and vanishes again at point B. So from this model, we can see 
that after removing the mean-field force, the effective short-range 
part of gravitational force of a particle is just dictated by its nearby 
particles. We separate the gravitational potential and the density as 
$(jc,t) = <bb(x,t) + S<b{x,t),a.ndL p{x,t) = pb{x,t) + Sp{x,t), 
in which ^b{x,t), pbix) are the long-range mean-field quantities, 
and 5<^{x,t), Sp{x,t) are the short-range residuals, respectively. 
With Poisson equation V^"l? — AnGp, these long- and short-range 
dynamical variables satisfy their respective Poisson equations 

V^$6 = inGpt, 

V^(5$ = AivGSp. (1) 

Thus, we can see that, in this mean-field model, the mean, 
smoothed-out gravitational potential is treated as a background 
field, corresponding to large-scale perturbations; whereas the short- 
range residual effects are regarded as fluctuations in this back- 
ground, corresponding to small-scale perturbations. In the mean- 
field model, a particle couples not directly with particles at large 
separations, but only with this background potential, with the resid- 
ual short-range effects determined by a few particles nearby. In 
later sections, we will see the great significance of this separation 
of gravitational force, as well as different physical implications by 
distinguishing between these two types of perturbations. 

2.2 Violent relaxation, two-body relaxation and Landau 
damping 

It is the instability of the mean gravitational potential, i.e. 
d^b /dt 7^ 0, that induces a strong, long-range relaxation mech- 
anism, called violent relaxation (.Lvnden-Be ll 1967). When violent 
relaxation ceases, i.e. d^b/dt = 0, it leaves behind for the system 
the virial equilibrium state, 2Ek / Ev = — 1, or equivalently, a hy- 
drostatic equilibriu m state, which is described by the Jeans equa- 
tion (see lBinnev & TremaindllOO^) 

dp 2/3p a* GMp 

:3~ ^ = ^p-^^ 2~' (2) 

dr r or 

where the 'pressure' p is defined as p = pa^, and is the radial 
velocity dispersion. /? is the velocity anisotropy parameter, usually 
defined as /? = 1 — aj^/a^, where is the one-dimensional tan- 
gential velocity dispersion, with = a]^ — a^. 

The short-range residual effects of the gravitational forces be- 
tween these particles are manifested in the following two forms of 
relaxation mechanisms. One is two-body relaxation, which arises 
from encounters of two bodies, and is only efficient locally, or at 
small scales of the system. At first, the system stays at its maxi- 
mum radius with vanishing kinetic (or thermal) energy, and then is 
driven by violent relaxation to approach the virial equilibrium state. 
It is well known that, at this stage, two-body relaxation operates in 



a far longer time-scale than violent relaxation, so that it should ac- 
tually be neglected when violent relaxation is functioning. 

The other short-range relaxation mechanism operating at 
small scales in a self-gravitating sy stem is gravi t ationa l Landau 
damping, which was firs t realiz ed bv lLvnden-Belll jl962h and fur- 
ther studied bv iKandrud ( Il998h . A brief introduction to this short- 
range relaxation is the following. A density perturbation with wave- 
length shorter than the Jeans length, Aj, will induce the Jeans 
waves; otherwise the perturbation will induce the Jeans instability. 
Such Jeans waves in self-gravitating systems, heavily damped by 
interactions with particles, are the counterparts of Langmuir waves 
in plasma. From the dispersion relati on of the self-gravitating sys- 
tems (equation 5.66 and fig. 5.2 of Binnev & Tremainal2008l) . we 
can see that Landau damping is also very efficient locally. Hence, 
the two short-range relaxation mechanisms are responsible for the 
local equilibria of the system at small scales. 

These relaxation mechanisms are picturized in Fig.fTJb). We 
can see that the three relaxation mechanisms, violent relaxation, 
two-body relaxation and Landau damping operate at different evo- 
lutionary stages and different spatial scales of the system: (1) two- 
body relaxation, which arises from interactions between particles 
themselves, is a short-range relaxation mechanism; (2) Landau 
damping, which is another short-range relaxation, arises from inter- 
actions between particles and the Jeans waves; (3) whereas if per- 
turbations take place in large-scale regions in the early evolutionary 
stage, such large-scale perturbations will induce long-range violent 
relaxation, i.e. the large-range violent changing of the gravitational 
potential. We will see later the significance of the discrimination 
between these three types of relaxation mechanisms. 

2.3 Entropy stationary state 

After violent relaxation terminates, the system stays in the hy- 
drostatic equilibrium state, described by the Jeans equation, equa- 
tion ([2}. In this section, we have interests to investigate whether 
there are some extra physics resulting from violent relaxation be- 
sides the hydrostatic equilibrium. 

As mentioned previously, there have been great disputes about 
the validity of the concept of entropy for self-gravitating systems 
for nearly five decades. Ignoring the controversy, we tentatively 
employ the entropy principle to explore the possible thermody- 
namic relations of self-gravitating dark matter haloes. 

We choose the specific entropy in the phenomenological form 
as the following 

s = ln(//^p-^''^), (3) 

which was first used bv lWhite & Narav"aril ( 1 19871) for the isotropic 
self-bounded gas sphere, but here we assume it can also apply to 
the velocity-anisotropic case. This expression is the standard en- 
tropy form of an ideal gas in terms of p and p, without the nor- 
malizing factor and an additive constant, and we borrow it for self- 
gravitating systems without any a priori reasons. With the specific 
entropy given, the total entropy of the system is 

St = 4-Kr^psdr = / 47vr^ p\n{p^^^ p~^^^)dr. (4) 
Jo Jo 

As usual, the constraints are the total mass and energy conser- 
vation. The mass function M(r) of a spherical gravitating system 
is 

M{r) = f 4nr'^p{r')dr', (5) 
Jo 
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and the total mass of the system is 
Mt = / 47rr^p(r)dr. 



(6) 



For the convenience of variational calculus, we set y{r) = M (r). 
So the mass-conservation constraint is converted to satisfying the 
following boundary conditions: 



y{r)\r=o = 0, y{r)\r=aa = Mt, 



(7) 



and we replace p{r) in the entropy functional with y'{r), the first 
derivative of ?/(r), as 



p{r) = 



1 dM(r) y'(r) 



(8) 



Airr^ dr 4:Tvr'^ 
The total binding energy of the system with spherical symmetry is 

J 2\r-r'\ 
= -47rG / M(r)p{r)rdr. 



(9) 



So by using the virial theorem, the total energy of the system is 



(10) 



The constrained entropy is 

St,c = St + \Et, (11) 

where A is the Lagrangian multiplier. Substitute equations lO, 
([8), and dlOl l into equation illl . we have 



St.c — 



F{r,y,y')Ar 



/, / 3/2/ y 



A 



Gyy' 
2r 



dr. (12) 



We treat p as a known function of r, and only calculate the first- 
order variation of 5*4, c with respect to y and y' . We obtain the fol- 
lowing Euler-Lagrangian equation from ^St^c = as 

^ = Alnfp3/2(^.-5/2 

dr dr V 'inr'^ 



dr 2r^ 



(13) 



Since this equation comes from the first-order variational calculus, 
in what follows, we call equation l |13t as the entropy stationary 
equation, or simply, the entropy equation. Additionally, from equa- 
tion it is interesting to notice that 



a — ln(p ' p 



3/2„-5/2n 



= -ln^. 



(14) 



That is, the so-called pseudo-phase-space density p/aj! seems to be 
related with the specific entropy s, which implies that the pseudo- 
phase-space density may re ally have a physical origin. This rela- 
tion has be en recognized by Isinnev & Tremaind ( l2008l see equa- 
tion 7. 121). IVass et al.l ( 120091) also suggested such a correlation be- 
tween them that the phase-space density profile likely reflects the 
distribution of entropy, which dark matter acquires as it is accreted 
on to a growing halo. 

With a simple dimensional analysis, we can write the La- 
grangian multiplier A in equation l ll3t as 



A : 



GpsTi 



(15) 



with p being a positive dimensionless number, so that A is always 
negative. In the following, we can see that equations ([2]l and ( 113b 
have no solutions if A is positive. We can always set /i = 1 by 
choosing appropriate values for the characteristic density ps and 
scale rs- 

In this work, we use the radial pseudo-phase-space density 
p/a'^. instead of the more general p/o'^ to formulate the specific 
entropy, which by no means indicates that the former is physically 
superior to the latter. In fact, we have tried many combinations of 
(7, but all failed to yield reasonable results. For example, take the 
form of a as, — arUea^ — a-'r{l -— /?), and perform the same 
variational calculus as for equation il3\ . we get the equation 



dln(p3/V/') 



d/3 



d In r 



-A 



GM 



(16) 



or with ct' 
following result 



(1 - P) dlnr 2r 

a'^{3 ~ 2/3)^/^ we derive the 



2 I 2 I 2\3/2 



dp 



-A 



GM 
~2r^' 



(17) 



dlnr (3 -2/3) dlnr 

Compared with the entropy stationary equation il3l . it can be 
seen that there appears an additional term, — d/3/dlnr(l — /3), or 
— 3d/3/dlnr(3 — 2/3) in the above two equations. From fig. 9 of 
Navarro et al. (2008), we can see that /? changes sharply at the out- 
skirts of the dark matter halo, so the absolute value of d/3/d In r at 
the outer boundary is huge. Such an additional term containing /3 
will greatly affect the behaviour of the solution of equations l ll6t or 
( I17t . so that we cannot even get a convergent solution. 

Moreover, there is also another problem that we cannot per- 
form the variational calculus with respect to the pressure p, since 
such a resultant equation is meaningless and unacceptable, so in the 
current work we have to treat p as a known function. As a conse- 
quence, such a treatment is inconsistent, and hence at this moment 
we can obtain only one equation by using the entropy principle. 
Because of this inconsistency, our results are incomplete. 

These problems indicate the great complexity of the statistical 
physics of self-gravitating systems, and imply that, either the en- 
tropy form may be even more complicated than the current forms 
we have already considered, or we may still lack some crucial 
physics, such as some extra macroscopic constraints besides the 
mass and energy conservation. We believe these problems should 
be the targets of future investigations. 



3 RESULTS 

Up to now, for the three independent variables p, p and /3, we just 
have two equations, the Jeans equation (|2]l, and the entropy equa- 
tion il3l . In order to obtain definite solutions of them, we have to 
provide an extra equation to close the equation system. This can be 
done by assigning one of the three variables with a known func- 
tion. In the following, we explore the solutions in two cases, re- 
spectively: (1) given the anisotropy parameter /3, and (2) given the 
density profile p. 



3.1 Given anisotropy parameter /3 

With the anisotropy parameter /? given, the density profiles can be 
solved, and these solutions are shown in Fig.|2] For simplicity, we 
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Figure 2. Density functions as tlie solution of equations (5) and )13t . with /3 
being a given function of r in equation U8t . Tlie fiducial solution of /3 = 
is indicated as 'Isothermal' in the figure; others are shown as /3 ^ 0. 



specify /3 as a step function of r as 



(18) 



in which C is arbitrarily specified as, say 0.5, and the values of 
the constant /?o are indicated in Fig. [2] First, we concentrate on the 
case of /? = 0. Here, with the Lagrangian multiplier A chosen to be 
negative, the two equations have neat analytic solutions as p oc 
and p (X r~^, which are exactly the isothermal solutions. It is well 
known that in such a solution, the system has infinite mass, energy 
and extent, which is similar to lLvnden-Belli l ll967h 's result, and so 
is not acceptable. 

Such an unconfine d density configurat ion of se l f-grav itating 
systems once misg uided lLvnden-Belll h96f) and Shu 1 1978) to be- 
lieve that violent relaxation is not a complete relaxation mecha- 
nism, and so i ntroduce the q uest for a theory of incomplete relax- 
ation (see also lMadserill 19871) . 

However, we can see that cases with /? > are different from 
the isothermal solution. These solutions gradually deviate down- 
wards from the isothermal solution and descend more and more 
rapidly with increasing radius. The larger the value of /?, the more 
decline of the solution, and the density functions drop down sharply 
at some outer radii. It is obvious to see that /? is of great signif- 
icance to make the solutions have finite mass, energy and spatial 
extent, and therefore, with /3 involved, it is not necessary to intro- 
duce a phase-space cut-off, o r incomplete relaxation mechanisms 
jLvnden-Bellll 1967l : IShij| 19781) . 



3.2 Given density profile p 

Now let us discuss the other case, i.e. solving equations (|2} and 
(|13 b bv specifying a dark halo de nsity profile p(r). H ere, we use 
the lEinastol j 19651) profile (also see iNavarro et alj200^ . 



ln{p{r)/ps 



'-[{r/r 



1], 



(19) 



with the s hape parameter a — 0.170, the same as the one taken by 
iNavarro e t al. (2008) for their Aq-A halo simulation series. 

We solve the two equations from the outer boundary to- 
wards the inner of the dark matter halo. The two solutions are 
shown as /3(r) and the radial pseudo-phase-space density p/crj? 
in Fig. [3] Compared with the corresponding simulation results in 
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Figure 3. (a) /3(r) as one of the solutions of equations and il3\ . with the 
density profile given as the Einasto profile in equation 1191 . We numerically 
integrate the two equations towards the inner of the dark matter halo, from 
the outer boundary at r = 16rs, where the logarithmic slope of the Einasto 
profile is dlnp/dlnr Ri —3.2. The initial values of /3 are shown in the 
figure. The arrow indicates the position of the local peak, at about r Ri 5rs. 
(b) The other solution of equations (2) and 1131 , shown as the radial pseudo- 
phase-space density profile p/cj? as functions of radius r. The vertical scale 
is unnormalized. 



iNavarro et al. I ( |2008| , see figs 9 and 12), we can see that our solu- 
tions in the inner part of the dark matter halo are completely un- 
acceptable. We ascribe this disagreement between our theoretical 
predictions and the simul ation results to the phenome non called 
gravitational degeneracy dCh avams &Sommerialll998l) . which is 
not included in this paper. We briefly introduce this effect in the 

follo wing. 

iLvnden-Belll l ll967h derived a distribution of self-gravitating 
stellar systems that is analogous to the Fermi-Dirac distribution. 
In his work, the gravitating particles are grouped into, and rep- 
resented by the coarse-grained, constant-density phase-space ele- 
ments, which have fixed and finite phase volume, and cannot over- 
lap to each other. It is this Fermi-Dirac-type distribution that leads 
to the existence of the gravitational degeneracy in the inner region 
of self-gravitating systems, where is much denser than the other 
parts. Such a gravitational degeneracy provides an exclusive effect 
against the gravitational coll a pse to avoid the gravothermal catas- 
trophe dLvnden-Bell & Wood|[l968l) . It should not be confused that 
the gravitational degeneracy results not from Pauli exclusion princi- 
ple of quantum mechanics, but from the conservation of the overall 
phase-space volume of these phase elements. In strong degener- 
ate region, i.e. the inner part of the dark halo, the close encounters 
should not be ignored so that the specific entropy, equation 
which is just suitable for ideal gas composed of point-like parti- 
cles, is not proper. It is valid only in the non-degenerate limit, i.e. 
the outer low-density region of the dark halo. 
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If the degenerate effect is included , the proper specific entrop y 
may be better generalized to, say, the White & NaravanI jl987h 's 



form 



s = -lnF-(-^-l)lri(7,-F), 



(20) 



where F = p/a^ is the pseudo-phase-space density, and rj is the 
maximum phase-space density, which embodies an exclusive ef- 
fect that prevents different phase elements from overlapping to each 
other in the final state. In the current work, however, we restrict our- 
selves not to include this complexity of gravitational degeneracy. 

According to the behaviour of the solutions, we draw a vertical 
dotted line at r = 2rs in Fig. [3] and an arrow pointing to the local 
peak at r ~ 5rs in the upper panel, to roughly separate the halo 
into the innermost strong-degenerate (with the radius r < 2rs), the 
middle weak-degenerate (2rs < r < 5rs) and the outmost non- 
degenerate (5rs < r < 16rs) region. As mentioned above, due to 
the ignoring of gravitational degeneracy, our theoretical predictions 
in the inner strong-degenerate region of the halo are not reliable at 
all. So in the following we just concentrate on the solutions in the 
outer regions. 

The velocity-anisotropy parameter for dark matter haloes that 
is greater than for baryon s has been observationally conformed 
jHansen & Piffarettil I200I iHost et al.l l2009h . whereas the latest 
high-resolution simulations revealed much more details that are 
absent in the observations. Nav arro et al ., (2008) and Lu dlow et al.l 
( l2010h found the non-monotonic radial dependence of /3, i.e. haloes 
are nearly isotropic near the center, radially biased to the maximum 
at some radius and approximately i sotropic again in the out skirts, 
which is against the earlier result of iHansen & Moor3 (I2OO6I. who 
found that the anisotropy parameter just increases monotonically 
outwards from the iimer part of the dark matter halo. 

From Fig. |3] we see that the predicted /3 agrees qua ntitatively 
well w ith the corresponding simulation results of Nav arro et al.l 
( l2008h in the non-degenerate region. In particular, it is worth em- 
phasizing that our theoretical predictions completely reproduce 
the non-monotonic radial dependence of /3 outside the strong- 
degenerate kernel. Besides, our predictions are also in qualita- 
tive agreement with the simulation results in the middle weak- 
degenerate region of the dark halo. 

Similar to the case of /3, the predicted pseudo-phase-space 
densities p/cr^ in the inner strong-degenerate kernel also greatly 
deviate from the simulated results, for the same reason as that for 
j3, i.e. ignoring of the gravitational degeneracy. In the outer regions 
with r > 2rs, however, the predictions, similar to the simulation 
results, can also be well fitted by p/ ff^ ~ ^.-1-875 _ jj^g power 
law predicted bv ISertschingej l ll985h 's similarity solution for in- 
fall onto a point mass in an otherwise unperturbed Einstein-de Sit- 
ter universe - except for an upward deviation in the halo outskirts. 
We are delight ed to see that our res ult is in good agreement with 
the finding by iLudlow et al.l | |2010|) . that such a curve-up devia- 
tion also exists in their simulations. This significant upward de- 
viation from the simple power law of p/af, ~ j-"^'^''^ indicates 
that the great complexity exists in equilibrated dark matter haloes 
dSchmidt. Hansen & MaccifilbOOSl) . 

Thus, we have seen that, compared with the simulations, 
our theoretical results provide good agreements in the outer non- 
degenerate region of the dark halo, and also acceptable agree- 
ments in the weak-degenerate region, respectively. We conclude 
that these agreements between our theoretical predictions and sim- 
ulation results sufficiently indicate that both the concept of entropy 
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Figure 4. Two modes of perturbation. The dashed lines in both panels indi- 
cate the equilibrium configurations, (a) Large-scale mass perturbation, i.e. 
A|(<5j;)^ 2> |-B|(i5j;')^; (b) small-scale density perturbation/fluctuation, 
i.e. \ A\ {&yf < \B\ (<5j/')^, where A and B are defined in equation WS . 



and the entropy principle should be still valid for long-range self- 
gravitating systems. 



4 GENERALIZATION OF THE SECOND LAW OF 
THERMODYNAMICS 

4.1 Second variation 

We would like to inspect whether the solution of the first-order vari- 
ational calculus, i.e. equation J13b . is a maximum, minimum or any 
other types of solutions. For this sake, we calculate the second- 
order variation. By straightforward substitution, the second-order 
variation of the constrained entropy of equation J12l ( is 



where 



A = 



[A{Syf + B{5y'f]Ar 



[A{SMY + B{A-Kr5py]Ar, 



(21) 



B 



cf£ _ d_f d'^F \ _ ^ G 
dy^ dr \ dydy' J 2r'^ ' 

d^F 



_5_ 

'2y' 



(22) 



For the convenience of the following analysis, we call the first term 
in the above integral, A{SM)^, as mass perturbation, and the sec- 
ond term, B{ATTr^Sp)^, as density perturbation. According to the 
result of the first-order variational calculus, we know A < 0, so we 
can immediately see that A > 0, and B < 0. This result indicates 
that the entropy stationary solution, equation l ll3t . is neither a max- 
imum nor a minimum, but a saddle-point solution of SSt,c ~ 0. So 
at first sight, it seems that the solution might be unstable, and the 
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entropy form of equation ^ may not be proper for self-gravitating 
systems. 

In fact, this apparent catastrophe contains a profound new 
physics that we may have never realized before. In what follows, 
we will analyze it in detail. 



4.2 The complementary second law of thermodynamics for 
self-gravitating systems 

First, we assume that 2> \B\(5y')'^, or equivalently, 

-X^{5Mf^^^{5p)\ (23) 

Here, the mass perturbation is predominant, whereas the density 
perturbation, i.e. the variation of the first derivative of mass func- 
tion, is negligible. As a result, such a mass perturbation must be 
large-scale, which is shown in Fig.|4la). One may ask what phys- 
ical processes can cause such kind of mass perturbation. A rough 
analysis shows that the dynamical processes, such as gravitational 
collapse or the system's radial large-range oscillations with large 
amplitude, can all lead to the mass perturbations, which must in- 
duce the violent changing of the long-range, background gravita- 
tional potential with respect to time, i.e. d^t/dt 0. This is ex- 
actly the feature of violent relaxation. 

So from equation | |2U . with A > and without small- 
scale density perturbation, we see that the entropy of the iso- 
lated self-gravitating systems, driven by violent relaxation, will al- 
ways automatically decrease to the minimum until violent relax- 
ation terminates. Furthermore, we notice that the mass perturbation, 
— AG((5A/)^/2r^, does not contain any localized information, say, 
some function of the density p(r), hence we can further identify 
this equilibrium entropy as the global minimum entropy. 

Our result is consistent with Binney's argument, which was 
designed to argue against the validity of Boltzmann-Gibbs entropy 
as well as the maximum entropy principle for self-gravitating sys- 
tems. This argument briefly states as follows. The total mass and 
total energy of the system are assumed to be Mt and Et (neg- 
ative), respectively. A large fraction of the mass and the energy, 
say AIt — p and Et — e, are concentrated in a bound static core, 
whereas the remaining mass p and thermal energy e are deposited 
in a diffuse outer envelope. As e — > the envelope becomes 
larger and larger in radius, but its density becomes more and more 
tenuous, such that the associated entropy will diverge to infinity 
teinnev & Tre maine 2008). So the fact that entropy is unbounded 
above is consistent with our finding^ 

Besides Binney's argument, lAntonovl ( [T962I) also proved that 
global maximum entropy state s do not exist for self-gravita ting sys- 
tems (see also appendix III of Lvnden-Bell & Woo j[l968 ). which 
is also a support to our conclusion. 

Now we turn to the analysis of density perturbation. Here, we 
have \A\(5y)^ < \B\{Sy')^, or equivalently, 

-X^iSMf«^-9^i5pf. (24) 
2r^ p 

In the context of density perturbation, the mass perturbation is neg- 
ligible, which implies d^t /9t — 0, i.e. there is no violent relax- 
ation. So we can see from Fig. |4{b) that the density perturbation 
is just small-scale perturbation. Such small-scale density perturba- 
tions will induce only short-range relaxations, i.e. dS^/dt 7^ 0, 
such as two-body relaxation and Landau damping. With B < 
and without mass perturbation, we can see from equation i2H that 



the entropy is the maximum for the system under short-range relax- 
ations. Moreover, because of the small-scale nature of the density 
perturbation, we can further identify this equilibrium entropy as the 
local maximum entropy, not for the whole entity, but just localized 
for every and any small part, or roughly, coarse-grain, of the sys- 
tem. 

We further illustrate the differences of the long- and short- 
range relaxation mechanisms in Fig. [5] When the system has expe- 
rienced a large-scale mass perturbation and departures greatly from 
the equilibrium state, long-range violent relaxation will be excited 
and drive the system back to the equilibrium state, i.e. the minimum 
entropy state. Nevertheless, if violent relaxation nearly or com- 
pletely terminates, then short-range two-body relaxation/Landau 
damping operate locally to maintain the local maximum entropy 
state for every/any small part of the system. In fact, in the usual 
thermodynamics with no long-range relaxation mechanisms, if the 
system is extremely large and unconfined, such as the atmosphere 
of the Earth, it is impossible to acquire a global, but just local ther- 
modynamic equilibria, for the system. So we can see that things are 
indeed complicated for the self-gravitating system, but our results 
do not conflict with the standard second law of thermodynamics at 
all, since the latter is just valid for the systems dominated by short- 
range relaxation mechanisms. We here just generalize the second 
law to accommodate the long-range, self-gravitating systems. 

To summarize, we re-express the second law of thermody- 
namics for self-gravitating systems in the seemingly paradoxical, 
but actually complementary way as follows: for an isolated self- 
gravitating system, if long-range relaxation (i.e. violent relaxation) 
is predominant, the entropy of the system never increases; whereas 
if short-range relaxations (i.e. two-body relaxation together with 
Landau damping) are predominant, the entropy never decreases. An 
additional explanation for this revision of the second law of thermo- 
dynamics states that the equilibrium states of self-gravitating sys- 
tems are the global minimum state for the whole system under long- 
range relaxation, but the local maximum entropy state for every and 
any small part of the system under short-range relaxations. It is only 
for the long-range self-gravitating systems that there exists such a 
minimum-maximum entropy duality, and the corresponding com- 
plementarity for the second law of thermodynamics. For the usual 
thermodynamic systems, there is no distinction between long-range 
and short-range relaxations, so there exists no such a complemen- 
tarity for the standard second law of thermodynamics. 

As addre s sed in the In troduction, some earlier authors (e.g. 
ISridhaJ 1 1 9871 : ISokeJ 1 19961) have already realized that the H- 
function, or entropy of self-gravitating systems, m ay not neces- 
sarily increase during the violent relaxation process. ISokeiHl996h 
even distinguished between two phases of violent relaxation, in 
which the ff-function may decrease with time in the first phase, 
i.e. the 'vigorous phase' they called, whereas in the second phase, 
i.e. the 'calm phase' they called, it is always a non-decreasing func- 
tion of time. If we tentatively regard Landau damping as the second 
phase of violent relaxation as they referred to, one can notice that 
their finding is very close to our result, which is a great support 
to the complementarity we found in this work. However, compared 
with our distinction between the three relaxation mechanisms, their 
two-phase classification seems slightly vague and inaccurate, and 
they did not relate these two phases to the long- and short-range 
relaxations, either. Most of all, it seems that they did not realize 
the new physics, i.e. the complementary second law of thermo- 
dynamics, underlying such a decreasing/incresing property of the 
_ff -function. 
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Figure S. Illustration of the effects of the three relaxation mechanisms. The horizontal dashed lines in both panels indicate the entropy in equilibrium states, (a) 
If the system deviates far beyond its equilibrium state, it will be driven by violent relaxation back to its equilibrium again, (b) Magnification of the local near- 
equilibrium state in panel (a). Under long-range violent relaxation, the system approaches its global minimum entropy state, whereas short-range two-body 
relaxation and Landau damping collectively operate to maintain this very entropy as the local maximum for every and any small part of the system. 
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Figure 6. A thought experiment helps to demonstrate the global minimum-local maximum entropy states of self-gravitating systems. Arrows indicate the time 
direction. From (I) to (III) are the free expansion of ideal gas; from (IV) to (VI) are the time reversal of the above procedure. 



We propose a thought experiment in Fig. |6] which may help 
to understand this generalized, or complementary, second law of 
thermodynamics. The processes are the following. (I) An ideal gas 
is confined within an adiabatic box, which resides in the vacuum. 
The short-range relaxations, namely, collisions of the particles be- 
tween themselves and with the walls of the box, maintain the ideal 



gas in its thermal equilibrium, i.e. the maximum entropy state. (II) 
Suddenly remove the walls, and the gas begins to expand outwards. 
(Ill) This is just the familiar free-expansion process of an ideal gas, 
in which the gas will expand all the way without ending, and we 
know that the entropy will increase to infinity along with the ex- 
pansion. Now we time reverse the above processes from (I) to (III), 



10 p. He and D.-B. Kang 



and then the time-reversal processes proceed in this way: (IV) the 
gas moves inwards to concentrate to a small spatial region, pulled 
by some fictitious 'long-range attractive force', and of course the 
entropy of the gas decreases with this process. (V) The ideal gas 
is suddenly enclosed by adiabatic walls, and then the 'long-range 
relaxation' stops, leaving behind a mechanical equilibrium of the 
ideal gas. The entropy of the system also stops decreasing, so we 
may say the system attains its minimum entropy state. (VI) The 
equilibrated ideal gas implemented by the 'long-range relaxation' 
is then maintained by the same short-range relaxations as in stage 
(I), so that the entropy is the maximum under the short-range re- 
laxations. We can see that, except for the sudden establishing of 
the potential well, i.e. the adiabatic walls, this time-reversal proce- 
dure greatly resembles the gravitational collapsing process under 
the action of gravitational force. This thought experiment helps to 
understand the seemingly strange behavior of the minimum - max- 
imum entropy states of self-gravitating systems. 



5 SUMMARY AND CONCLUSIONS 

The statistical mechanics of self-gravitating systems has long been 
believed as an inviable discipline. The controversy mostly focuses 
on whether the concept of entropy is still valid for long-range inter- 
action systems, such as self-gravitating systems. In this work, we 
perform a preliminary study to this age-old problem. We summa- 
rize our findings in the following. 

When self-gravitating systems attain their equilibrium states, 
the systems generally possess some geometrical symmetry, such as 
a spherical shape. Under this condition, it is desirable to approx- 
imate the self-gravitating systems by the mean-field model. That 
is, the long-range part of gravitational forces between the parti- 
cles is collectively manifested as a smooth background force field, 
with the short-range residual effects determined by nearby parti- 
cles. Such a separation of long- and short-range gravitational force 
breaks the long-range coupling of a particle with the others at all 
different range of scales, and definitely renders the problem more 
manageable. 

Entropy is a state function that is introduced to describe the 
degree of order or chaos of a thermodynamic system. This concept, 
however, is of great significance in statistical mechanics, not just 
because it is a thermodynamic state function, but also because it 
constitutes a variational principle called the entropy principle (with 
the obsolete name 'principle of maximum entropy'), from which 
the equations or relations of thermodynamic quantities can be de- 
rived. It plays the same role for statistical mechanics as the one 
played by the principle of least action for fundamental dynamics. 
Without an appropriate definition of entropy of the thermodynamic 
system under consideration, it is impossible to establish the statisti- 
cal mechanics for this system, as has been exactly the case of self- 
gravitating systems for the past several decades. In this work, we 
tentatively employ the phenomenological entropy form of ideal gas, 
first used by White & Narayan ( 1987), to explore the possible ther- 
modynamic relations of the equilibrium states of self-gravitating 
systems. 

By calculating the first-order variation with respect to the den- 
sity function p, subject to the constraints of mass and energy con- 
servation, we derive an entropy stationary equation. However, we 
acknowledge that our treatment is somewhat deficient in complete- 
ness or consistency, in that we have to treat the pressure p as a 
known function, and cannot calculate the variation with respect to 
it. As a result, we can derive only one equation, i.e. the entropy sta- 



tionary equation. Mo reover, we have also tried some other entropy 
forms other than the IWhite & Naravaij fT987h form, but all failed 
to yield any reasonable results. These problems imply that either 
the entropy form may be even more complicated than the current 
forms we have already considered, or we may still lack some cru- 
cial physics, such as some extra macroscopic constraints besides 
the mass and energy conservation. However, despite this incom- 
pleteness/inconsistency, when incorporated with the Jeans equa- 
tion, we do find some interesting results with this entropy equation. 

In the case of given anisotropy parameter /3, we solve the en- 
tropy stationary equation and the Jeans equation to obtain p and p. 
If f} — 0, the two equations possess the isothermal solutions. How- 
ever, solutions with /? > deviate downwards from the isother- 
mal solution, descend more and more rapidly with increasing ra- 
dius, and drop down sharply at some outer radii. We conclude 
that j3 is of great significance to make the solutions have finite 
mass, energy and spatial extent, and hence it is not necessary to 
introduce a phase-space cut-off, or incomplete relaxation mecha- 
nisms. According to the behaviour of the solutions, we separate 
the halo into the inner strong-degenerate (r < 2rs), the middle 
weak-degenerate (2rs < r < 5rs) and the outer non-degenerate 
(5rs < r < 16rs) region. In the case of given density profile p, say 
the Einasto profile, we derive the solutions of /3 an d p/cr^, which 
are in good agreements with the simulation results of lNavarro et al.l 
(2008) in the outer non-degenerate region of the dark halo, and 
also in rough agreements with the simulations in the middle weak- 
degenerate region. In particular, it is worth emphasizing that our 
theory predicts an upward deviation of p/a^ from the power law, 
p/crj! ~ ).. -i.875^ ^^^^ ou tskirts, which is consistent with the 
finding bv lLudlowetai] ( l20I0h . that such a curve-up deviation also 
exists in their simulations. However, we cannot explain the simu- 
lation results in the inner strong-degenerate region, since we have 
not considered the effect of gravitational degeneracy in the current 
theory. 

By calculating the second-order variation of the entropy func- 
tional, we find that the stationary solution of the first-order varia- 
tion is neither a maximum nor a minimum, but a saddle-point so- 
lution. We argue that such a saddle-point solution by no means in- 
dicates a failure of the theory, but rather suggests a new physics 
that we may have never realized before. Inspired by this saddle- 
point solution, we distinguish between two types of perturbations 
in self-gravitating systems, namely, the large-scale mass perturba- 
tion and the small-scale density perturbation, which correspond 
to long-range violent relaxation and short-range two-body relax- 
ation/Landau damping, respectively. These two types of relaxation 
mechanisms operate in different fashions. For an isolated self- 
gravitating system, the entropy can automatically decrease under 
long-range violent relaxation, until it terminates and the system 
attains the equilibrium state with the global minimum entropy. 
This result is consistent with Antonov's proof, or Binney's argu- 
ment that there are no global maximum entropy states for self- 
gravitating systems. However, for every and any small part of the 
system in its equilibrium state, the entropy is actually the local max- 
imum under short-range relaxations, which is not in conflict at all 
with the standard second law of thermodynamics. This minimum- 
maximum entropy duality constitutes the complementary second 
law o f thermodynamics for long-range self-gravitating systems. 
ISokerl ( 1996) found the decreasing/increasing property of the gen- 
eral iif -functions in the two-phased violent relaxation, which can 
be regarded as a support to our finding in this work. 

Although we have not proposed a complete entropy form, it 
is extremely intriguing that this simple, or somewhat crude, en- 
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tropy form of ideal gas, by IWhite~& NaravanI ( Il987l) , can lead to 
such a profound new physics. Nevertheless, our study has defi- 
nitely demonstrated that the concept of entropy is still valid for 
self-gravitating systems. However, the standard second law of ther- 
modynamics should be revised, and the term 'principle of maxi- 
mum entropy' should be abandoned, and should be simply replaced 
by 'entropy principle'. The validity of entropy for self-gravitating 
systems means that the statistical mechanics and the relevant ther- 
modynamic relations can also, in principle, be established for such 
a long-range interaction system. We believe that our findings, es- 
pecially the complementary second law of thermodynamics, are 
interesting, and may offer important clues to the development of 
statistical mechanics for self-gravitating systems as well as other 
long-range interaction systems. 
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